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The degenerate Blume-Emery-Grifliths model for martensitic transformations is extended by in- 

- - - eluding both structural and magnetic degrees of freedom in order to elucidate premartensitic effects. 

O* , Special attention is paied to the effect of the magnetoelastic coupling in Ni2MnGa. The micro- 

scopic model is constructed and justified based on the analysis of the experimentally observed strain 
variables and precursor phenomena. The description includes the (local) tetragonal distortion, the 
amplitude of the plane-modulating strain, and the magnetization. The model is solved by means of 
S ' mean-field theory and Monte Carlo simulations. This last technique reveals the crucial importance 

^ I of fluctuations in pretransitional effects. The results show that a variety of premartensitic effects 

may appear due to the magnetoelastic coupling. In the mean-field formulation this coupling is 
quadratic in both the modulation amplitude and the magnetization. For large values of the magne- 
toelastic coupling parameter we find a premartensitic first-order transition line ending in a critical 
point. This critical point is responsible for the existence of large premartensitic fluctuations which 
manifest as broad peaks in the speciflc heat, not always associated with a true phase transition. 
The main conclusion is that premartensitic effects result from the interplay between the softness of 

Cq ' the anomalous phonon driving the modulation and the magnetoelastic coupling. In particular, the 

premartensitic transition occurs when such coupling is strong enough to freeze the involved mode 

\^ , phonon. The implication of the results in relation to the available experimental data is discussed. 
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I. INTRODUCTION 

Many metals and alloys undergo a so called martensitic transition (MT) from an open cubic phase at high tem- 
peratures to a more closed packed phase at lower temperatures Q. It is a displacive, diffusionless, first-order phase 
transition, accompanied by incomplete softening of certain transverse phonon modes. For Zr, which belongs to an 
ideally simple class of martensitic materials, the pure group IV-metals [^P, it was demonstrated that the first-order 
character can be understood as an effect of a coupling between two simultaneous strains: an internal two-plane shuffle 
strain and a uniform strain 0] . 

A rich variety of precursor phenomena ||5|] have been observed in (weakly) first-order MT. Some of them, as the 
intermediate tweed structures MM, are not common to all materials but others, intimately related to the transition 
mechanism, are present in almost all bcc systems studied so far. The most significant is the anomalously low Tyl2[110] 
phonon branch {[110] propagation, [110] polarization]), accompanied by a low value of the elastic constant C = 
(Cii — Ci2)/2. Moreover, both the phonon branch and the corresponding elastic stiffness soften with temperature. 
Recently, a lot of interest has been focused on the intermetallic Ni-Mn-Ga alloy close to the stoichiometric composition 
Ni2MnGa. It is the only known ferromagnetic fee Heusler alloy undergoing a MT on cooling. Besides its theoretical 
interest, it may be of technological importance too since it opens the possibility of controlling its shape memory 
properties (intimately related to the MT) by applying an external magnetic field [p[. It has a transition from a 
bcc (neglecting the atomic order) to a low temperature tetragonal phase bet, which is modulated by a five-plane 
shuffle strain |P,M. Particularly intriguing is that for nearly stoichiometric composition the full MT is preceded 
by an intermediate phase in which apparently only the shuffle strain is activated, but not the tetragonal strain 
[0,12 13 14. This intermediate phase consists in a micromodulated domain structure, without resulting macroscopic 



tetragonal deformation so that the cubic symmetry is preserved |14| . This is accompanied by a significant, although 
not complete, softening of the TA2 phonon branch at a wave vector ^0 = 0.33. Only at lower temperatures, at the 
martensitic transition point, the homogeneous tetragonal strain is activated (and the modulation changes slightly). 
This particular behavior observed in Ni2MnGa seems to be related to the influence of magnetism. The point is 
whether or not this intermediate modulated structure is a genuine new phase. Different behaviors have been reported 
in the literature. For some samples [ pl] , p"3yi^ ] there exists evidence for a true phase transition of very weak first-order 
which is driven by a magnetoelastic coupling. The main proof of that is the fact that the intermediate transition (IT) 
shifts with the external applied field [|ll|] w hile no (significant) magnetic field dependence has been found for the MT 
temperature ]|lq,tt6[. In other studies [|17|], the authors could not find any indication for an IT although precursors, 
clearly related to the magnetization of the sample, have been observed. Apparently, the only relevant difference in 
the samples used by Planes et al [O], Zheludev et al pqjlj] and Stuhr et al [O is the content in Mn. Very recently, it 
has been suggested [|| that the tetragonal phase can be suppressed by increasing the concentration of Ni at expenses 
of the content of Mn. Moreover, the magnetoelastic effects have been confirmed from the experimentally observed 
dependence of the elastic constants on an external magnetic field ]|l^ . In spite of the experimental evidence for the 
magnetoelastic coupling in Ni2MnGa, its microscopic origin has not been established yet. 

In the present paper we have theoretically investigated the nature of the bcc to bet transition and constructed 
a model in order to solve the puzzling behavior and to elucidate the role of the magnetic coupling, which counter- 
intuitively seems to couple the ferromagnetic order stronger to the modulating strain (77) than to the uniform tetragonal 
strain (e). We shall demonstrate that the situation can be described by the degenerate Blume- Emery-Griffiths model 
(DBEG) ]|l^ extended to include coupling to magnetic degrees of freedom, and with a interpretation of the variables, 
appropriate to the present case. 

The plan of the paper is the following: First, in section II, we provide an analysis of the experimental facts and a 
theoretical explanation of the observed phenomena, which is used to justify a microscopic model presented in section 
III. The model is studied using first mean-field theory (section IV) and next Monte Carlo simulation (section V) . The 
discussion from the comparative study is presented in section VI where we provide also our conclusions taking into 
account the available experimental data. 

II. EXPERIMENTAL FACTS AND THEORETICAL INTERPRETATION 

A. Experimental Facts 



The structural properties of Ni2MnGa have been investigated in a series of papers [Pl Jl 0| Jl l| , |l 2| , |1 3| , |l 4|l 7y20|] . At 
high temperatures the alloy has the fee (i2i) Heusler structure which, neglecting the atomic order, can be regarded 
as a bcc lattice. It is paramagnetic at high temperatures with the magnetic moments mainly on the Mn sites 
pl|. At temperatures below T^, it orders ferromagnetically with no particular easy direction of the moment. At 



the temperature Tm (< Tm) there is a (first-order) structural phase transition of the martensitic type to an average 
tetragonal structure, which additionally is modulated by a transverse five-layer shuffling strain. Prior to this transition 
precursor structures of tha.t phase as well as of the fee having a six-plane modulation have been observed in neutron 
scattering experiments [|l4| , [f7[ . This may happen also as a transition (first-order) at a temperature T/ {Tm < Tj < T„i) 
giving rise to a genuine intermediate phase |23| without any macroscopic tetragonal deformation urn . 

The above temperatures T,„, Tm and T/ are extremely sensitive to the composition and atomic ordering of the 
sample. Thus Tm may vary between 360K to 395K, whereas Tm may vary from 175K to 450K. In the sample studied 
by Stuhr et al. [|7J, Tm ~ 364K and Tm ~ 284K with precursor signs for T > Tm, but no intermediate phase was 
found. In the sample studied by Zheludev et al. [|l3|J l4f| Tm ~ 380K and Tm ~ 220K and an intermediate phase is 
observed below Tj « 260K. Similarly, Planes et al |11|] studied a sample with Tm ~ 381K and Tm ~ 175K and found an 
intermediate transition at Tj « 230K. In spite of the large variation in the temperatures, we assume the basic physics 
is the same, and hence we shall use the information obtained by observations on different samples. In particular, we 
shall attempt to explain the precursor and the intermediate phase phenomena by constructing an effective microscopic 
Hamiltonian, which allows analysis beyond a mean- field or Landau expansion treatments P, pl] , p3[| . 

B. The two-strain model 

Although Ni2MnGa is a metallic alloy for which both the structure and the magnetism is determined by the 
conduction electrons, it is instructive to consider a model system for the mismatch between cubic crystals, similarly 
as it is found in KCl grown on NaCl (001) ||2J,^. In this case, the very large (~ 17%) lattice constant mismatch 
makes the interface to buckle simultaneously in the [110] and [110] directions. This produces superstructures consisting 
of seven NaCl and six KBr layers, or multipla thereof, perpendicular to the interface, while preserving the square 
symmetry. Both the NaCl and the KBr crystals are modulated since they have similar elastic properties. This 
gives rise to superstructure peaks in the X-ray and Helium scattering spectra, precisely as those observed in the 
intermediate phase of Ni2MnGa, where one also would expect more higher order satellites in off-symmetry directions 
to be found, but not so far looked for experimentally. The buckling gives rise to a variation of the local [001] direction 
and a change in lattice constant perpendicular to the interface (corresponding to a local tetragonal distortion). We 
stress that for the ionic crystals the forces by no means triggers the local modulating strains, which are caused by 
the forced contact between unequal crystals. For Ni2MnGa the situation is rather different. Here, a nesting feature 
of the Fermi surface causes a strong electron-phonon coupling and an incipient soft phonon mode at q = (^,^,0) 
positions in the fee phase, where ^ ~ -^ ]|l4| , [17| . This is presumably the driving mechanism in Ni2MnGa, and giving, 
as a consequence, the tetragonal distortion. As precursor phenomena, quasi elastic peaks are observed at ^ = i and 



^ = i 1 17 1 corresponding to a similar six-plane modulation as discussed for the alkali salt interface. To match this to 
the fee (001) lattice plane it is advantageous to make a lattice mismatch and expand the lattice and to rumple the 
interface, or in other words making the [GDI] direction fluctuate in epitaxial-like angles. Stuhr et al [ p^ have shown 
that the six-plane fee modulation corresponds precisely to a five-plane modulation of the tetragonal phase, and the 
latter is found as a precursor phenomenon. An 5:6 expansion would be very drastic. However it suffices to create 
a superstructure of a common divisor, for example a 30 or 60 plane repeat distance. The latter would correspond 
to a lattice mismatch of 1.67%, which is very close to the mismatch observed between the fee and tetragonal phase 
of Ni2MnGa: ~ 1.6% [£1[. Hence we argue that the electronically driven six-plane modulation in turn also causes 
the tetragonal distortion. The theory for why it is advantageous for mismatched crystals at epitaxial interfaces to 
develop mutual superstructure peaks was discussed in more detail by Vives and Lindgard [ p6[ . The modulation occurs 
simultaneously in the q= [-^,^,0] and q= [i,— ^,0] directions, thus preserving the square symmetry of the (001) 
plane and yielding a modulated (001) plane, rather than a direction. In the fee structure there are three such equivalent 
planes (100), (010) and (001). In order to (essentially) preserve the volume, V, a change of which is expensive for 
the electrons, the c-axis [001] must shrink (from T = 4.2 K to 295 K it is indeed observed that V^ only increases by 
0.5% |2|]). Therefore a model must include a coupling between the plane modulation rj and the tetragonal strain e. 
Then, the question is: how can there be apparent separate temperatures for the onset of ordering of the two kind of 
strains? This is possible if in the intermediate phase the tetragonal strain is only local, varying in direction and also 
in the allowed directions in the fee crystal. Moreover, only in the martensitic phase microcrystals with a resulting 
tetragonal strain should be formed. The corresponding tetragonal structure is observed as highly mosaic with a large 
variation of the [001] directions [|l^], in agreement with the above picture. 



C. The Landau models 

Very recently, Landau models for the MT in cubic ferromagnetic materials have been proposed [ppq] ■ In these mod- 
els, the magnetoelastic coupling between the uniform strain tensor and the (vector) magnetization is fully considered. 
Nevertheless, they seem more appropriated to the study of the magnetic properties of the martensitic phase rather 
than to the analysis of the IT itself. Here, accordingly with the discussion given above, we shall adopt a different 
strategy and study the structural transition in Ni2MnGa in terms of a Landau expansion of the only most relevant 
strain and magnetization variables, including non-uniform strains or modulations. Then, similarly to the case for the 
bcc to hep transition in Zr Q, we include the following terms, which are allowed by symmetry: 

+ ^e^f, (1) 

where 77 is the discussed plane modulation strain and e the local uniform contraction perpendicular to that plane, but 
we do not consider higher order uniform strain terms. Here w^ — a{T — Tg) is the squared frequency |2^ for the 
incipient soft mode phonon with q = (i ^,0), other constants are positive, and C" = hiCn — C12), which is small 
and temperature dependent. By eliminating the local tetragonal strain e we can [Q, write the free-energy along the 
optimum energy path involving both e = — A^/ C" i]^ and r] as 

B = B-2A^/C'. (2) 

The coupling between the two strains therefore makes it possible for B to become negative and hence to cause a 
first-order transition at Tm before the soft mode transition at T = T^, even without the coupling with the magnetism. 
Next let us consider the influence of the magnetism. A ferromagnetic moment m can influence an itinerant magnet 
as Ni2MnGa in two ways. Firstly, giving rise to a splitting of the electronic energy bands, which is proportional to 
the amplitude \m\, but which is not sensitive to the moment direction. This will cause a change in the soft phonon 
frequency below T^, proportional to |mp: 

cj2(m) = a{T - Ts) - u\m\^ = a{T - T,) + biT - Tm), (3) 

where a, h and u are positive constants, and we have assumed a mean-field behavior for |mp. A measurement of the 
soft mode frequency-squared should therefore exhibit a kink at the magnetic ordering temperature. This is exactly 
the behavior observed by Stuhr et al [0, who found a — 0.018 meV^/K and b = 0.020 meV^/K. Unfortunately for 
the samples showing the IT |0,0,Q measurements in the paramagnetic phase are not available. Nevertheless, in 
what follows, we shall assume the behavior discussed above for all the samples. Hence the magnetic free-energy part 
can be written as: 

•^mag(m, ry) = - a Imp + - /3 |m|4 + • • • + ^ 7 \m\''f, (4) 

where a = A(Tc — T) and Tc is the magnetic transition in absence of magnetoelastic coupling. A and j3 are positive 
parameters and 7 = — /iu is yielding the coupling between the amplitude-squared of the magnetization and the effective 
modulating strain. By eliminating \m\'^ one can write an effective free-energy in the form of the eq.fl). This yields a 
further temperature dependence oi ujI and B, as discussed by Planes et al ||ll|] . 

The other possible coupling between the magnetization and the structure is by magnetostriction, which deforms 
the crystal in the direction of the magnetization ||8| . Experimentally, the easy direction of magnetization is not known 
with certainty, even in the tetragonal phase. Webster et al |£1[ propose it might be in the (111) directions of the L2i 
phase |2^ , but that other directions are almost as likely. Analysing their data perhaps allows the conclusion that it is 
at least not in the tetragonal [001] direction, which would normally have been the obvious choice. If so, it is confined 
to the (001) plane, say, which has four-fold symmetry and therefore not yielding a strong easy axis. The lowest 
order coupling would then be of the form -\-fh?e^ . An interesting possibility is if the easy direction is along the [100] 
direction, because this may belong to two different modulation planes, (001) and (010), therefore yielding a minimum 
magnetostriction because it cannot distinguish between a [001] and a [010] tetragonal strain. The effect is further 
reduced because of the equivalence between moments in the [100] and [010] directions in the modulated plane (001). 



Inclusion of the coupling can be done in the eq. (Q) , with no qualitative change (except for an induced dependence on 
the magnetization of C" ||ll|]), and hence we shall neglect it in the following. If the coupling is sufficiently strong, it 
would no doubt prevent the existence of an intermediate phase in Ni2MnGa. 

The conclusion of the above analysis is that the phases may be characterized by the minimum path modulation 
strain fj in eq.(p|) which includes a finite, but local tetragonal strain which is on average zero because of fluctuating 
directions. This is consistent with the observation of a significant broadening of the fee (002) peak in the neutron 
scattering data p7[ |. This strain is coupled to the magnetization such that a ferromagnetic moment on neighboring 
sites will favor a modulation also at these sites, due to a change in the local band structure, but without differentiating 
between the three possible planar modulations and accompanying tetragonal local strains. In order to give meaning 
to the notion of a local band structure, we consider sites not as the Mn atom positions, but at a more coarse grained 
level, e.g. at the unit cell level, say. At the MT a resulting tetragonal strain, e, may arise as the interfaces between 
the modulated fee and the modulated and locally tetragonal crystallites grow larger and thicker. This strain is on the 
average in the (001) directions but with considerable variations in the epitaxial angles relative to those. To construct 
a simplified model, we shall consider one with only two variables such that each has only two degrees of freedom. 
One variable with only two values is representing the plane modulations 77, and the other also with only two values 
representing the resulting tetragonal strain e. This does not correspond to considering the behavior of a subspace of 
variants, but constitutes a simplified statistical analogue to the physics of Ni2MnGa. 

III. THE MODEL 

The desired degrees of freedom can be represented by the p-degenerate Blume- Emery-Griffiths Hamiltonian (DBEG) 
|l9| . It was first introduced with the aim to account for the entropy stabilization of the high temperature phase in 
martensitic transitions. Recently, it has been shown |g^ that it is equivalent, with respect to universal properties, to 
that of the ordinary BEG model with the crystal field shifted by a term fc^rinp. 

Although the DBEG model is very simple it includes most of the relevant physical ingredients to understand MT, 
namely a multivariant low temperature deformed phase and a high temperature average cubic phase with enhanced 
entropy. The model is an extension of the ordinary three-state Blume-Emery-Griffiths Hamiltonian defined on a 
lattice, which we shall take as simple cubic (or square), as motivated above. On each lattice site i = 1, ..,N a. variable 
tJi = 1, 0, —1 represents the deformation state near each site on the lattice. The state a.i = represents the undistorted 
phase, and it is chosen to be p-fold degenerate (p > 1), in order to approximately account for the high entropy of 
vibration of the cubic phase. The states Ui = ±1 represent the distorted phase. The Hamiltonian accounting for the 
energy gain in having the same structure on neighboring sites was written as pSJ ] 

n.n. n.n. 

Km = -JY.'^^'^i - ^E(1 - '^»')(1 - ^^y (5) 

(«.i> {id) 

where the sums are performed over all nearest-neighbor pairs. In what follows we will take J > as the unit of energy 
and work in terms of reduced magnitudes defined as H* — "H/J, K* — K/J .... 

Two order parameters can be defined, 0i = X) ^i/^ and (j)2 = ^ ^f f^- The model in (||) was solved, for K* > 0, 
by mean- field and Monte Carlo simulation techniques |l3|. We found a phase transition from a cubic (disordered) 
phase with (j)i — to a tetragonal (ordered) phase with (f>i ^ 0. The behavior of the secondary order parameter 
02 in the disordered phase depends on K* and T*, but exhibits non-analytical behavior at the same temperature as 
(pi. Large values of K* stabilize the cubic phase. Moreover, the K* parameter and the degeneration p of the 0-state 
control the order of the transition, which changes from being of second order (for low values of K*) to first-order. 



The DBEG model |19| was introduced as a simplification and a further generalization of the Lindgard-Mouritsen 
model H], initially designed to mimic the bee to hep transition in Zr. Consequently, it naturally inherits the identifi- 
cation of the two order parameters suitable for Zr: (f>i is the two-plane shuffle strain and (J32 the homogeneous strain. 
In the previous work |19[ this was not emphasized since both <j>i and (j>2 exhibit a phase transition at the same point. 

For Ni2MnGa it is more natural to identify the order parameters in the opposite way. Therefore, in the present 
work, we define: 



The breaking symmetry order parameter e corresponds to the tetragonal distorsion in the sense that if e = (equal 
population of at = +1 and —1) it corresponds to having all the variants equally populated, and hence an average 
cubic phase. For q the relation is more involved because of the complicated physics of the modulating strain ij. In 
the high temperature cubic phase, when the ct^ variables distribute at random q — qq = 2/{p + 2). Let us assign 
the difference qo — q with the amplitude fj of the plane modulating strain without distinguishing between the three 
possible modulating planes. 

We now include the magnetic degrees of freedom by means of spin variables Si — ±1 (defined on the lattice site 
i = 1, --jN) having a ferromagnetic Ising interaction. Thus, the purely magnetic contribution is: 






- J* 



E 



Ji Oj 



(8) 



where J^ > 0. 

The total Hamiltonian should further include a coupling term between the structural and magnetic variables. We 
have argued in section II that the magnetic influence of electronic properties gives rise to a coupling between the 
magnetic moment and the plane modulation. On a microscopic level let us assume that the presence of a moment 
of neighboring sites gives rise to a modulation on neighboring sites as motivated previously. To describe this, let us 
consider the following symmetry allowed |pO| interaction contributions: 
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(id) 
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that may be rewritten as: 
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In the pure ferromagnetic phase {Si = 1) this Hamiltonian can be viewed as an Ising Model for the variables 1/2 — af. 
Thus a phase transition exists for {Uii + Uqq — 2Uiq) > and Uii = C/qo- ^'^^ simplicity, in what follows we shall take 
Uii = Uqq = and denote U* — Uiq < 0. Then, the coupling Hamiltonian becomes: 



HL,^2U*l^ S,S, (i-a^ 



^ (id) 



(11) 



As we shall see in the next section, in the mean-field approximation, the first term becomes of the form of the coupling 
term in eq.(H) whereas the last term gives a simple modification of the purely magnetic interaction J*^ defined in 
eq.(|2). Furthermore, ( |ll[) shows that it may be particularly convenient to choose p — 2, which gives qo — 1/2. 
The total Hamiltonian model for Ni2MnGa can then be written as 



n* = n 



M 









(12) 



with H^j, Ti'^ and Ti*„f respectively given by expressions (|5|), (||) and (pT|). We shall demonstrate that it is possible 
to split up the structural transition into one determined by the order parameter q, which we will associate with the 
IT and another one, determined by e, to be associated with the tetragonal deformation occurring at the MT. 



IV. MEAN-FIELD TREATMENT 



In this section we solve the presented model ( |l2| ) by using standard mean-field techniques. The state of order 
of the system depends on the occupation numbers N^ . This stands for the number of points in the structural 
a = —1,-1-1,0 and magnetic S = +,— state. There are six different occupation numbers which should fulfil the 
following normalization condition: 



N+ + N^ + N+ + N^ + N+^ + N: 



N, 



where N is the total number of points in the lattice. We define the following order parameters: 

i 
i 

Nm = Y^ S, = N+ + N+^ + N+ - TVf - NZ^ - N^ 

i 

Nmo = ^(1 - a,^)S, - N+ - N^ 

i 
i 

The corresponding entropy can be written as: 



SMp/kB = In 



Nl 



—,P 



<.N++N-) 



where p > 1 is the degeneracy factor of the 0-state. The mean-field expression of the free-energy per particle is: 

^liF ^ = -{e'+K*{l- qf + J*^m^ + 2t/*mo(m - mo)) 



(13) 
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- q- mo 



(20) 



where T* = ksT / zJ , and z is the coordination number of a given site. Notice that mi appears only in the entropic 
contribution to the free-energy. Standard minimization with respect wi renders the following relationship mi = 
e(m-mo) ^ which has to be fulfilled at all temperatures. Then, after substitution in eq.(20) we obtain the following 
expression for the free-energy as function of e, g, m, mo'. 

HiF = -(e^ +K*{1- q)^ + J:„m^ + 2U*mo{m - mo)) 

+ T*[(, + e)ln(i±l)+(,-.)ln(i^ 

, . , /g-|-m-mo\ f q - m + mo 

+ [q + m, — mo) m I + [q ~ m + mo) In 



+ (1 - q + mo) In + (1 - 9 - mo) In 



2(l-q)ln(|)-2glng 



(21) 



Further minimization with respect the other four order parameters yields to the next set of coupled equations 

^lni±i 
2 q~e 



-?lni±i, (22) 



^P^ 



'2' 



T* f {q + e){q - e){q + m - mQ){q - m + mo) 



K*il~q)+T*Hi^) =Vln r^ ,.......,„..,..... „ , „.,. . ^23) 



2 V g2(-J^ _g_|_„jjj)(-2 _g_^„p) 



J:.m + U*mo = ^ln '^^^-^\ (24) 

2 q — m + mo 

C/ (m - 2mo) = — In -^— — -. (25) 

2 (g + m — mo)(l — 9 — mo) 

Their solution gives the temperature dependence of the order parameters. Between all possible solutions only the 
absolute minima correspond to thermodynamic equilibrium. This requires the analysis of the second derivatives of 
the free-energy. 

The space of the model parameters of interest here is limited by the conditions Jj^ > 0, C/* < and J*^ + [/* > 0. 
Then, for appropriated values of K* ^ there exists three phase transitions at the temperatures T^^ < Tj < T*-^ 
associated with e, q and m. respectively. In what follows we shall fix the value of the magnetic interaction to Jj^ — 4.0, 
thus determining the distance between T^.j and T*^, and use different values of the coupling parameter U* {0 < U* < 
—J^ = —4). Then, the values of K* for which the IT exists are determined by eq. (^3|). Indeed, by setting e = 0, 
it follows that for K*{T*) = —2T*ln{p/2) the order parameter q has a continuous phase transition. The results will 
be presented for two values of the degeneracy factor, p = 2 and p — 4. In Table I we give the identification of the 
different phases in relation to the problem of interest here and their corresponding abbreviated notation. Figures 
n and show the temperature behavior of the order parameters for a given value of U* = —3.50 along the path 
determined by the condition K*{T*) = —2T* ln(p/2) (coexistence line). In particular, for p = 2 one has K*{T*) = 0. 
For the sake of clarity, mo(T*) is not shown. At high temperature, a magnetic transition appears at T^ — J^ + U* /2 
from a paramagnetic cubic phase (PC) (e = 0, g = 1/2, m ~ mo = 0) to a ferromagnetic pure cubic phase (FC) 



(e = 0, (7 = 1/2, m, = 2mo 7^ 0). From equation (23) it is easy to see that, for q = 1/2 (T* < T/), m, and mo are not 
independent but m, = 2mQ. At lower temperatures, the order parameter q separates into two branches q~^{q > 1/2) 
and q^{q < 1/2). This occurs at the critical point T/^ and separates two different ferromagnetic phases both with 
e = but with different values of q: a pure cubic phase (e = 0,^ = 1/2, m = 2mo 7^ 0) and an average cubic 
phase {e — 0,q ^ 1/2, m ^ mo 7^ 0). The two branches of q are identified with the tetragonal-like modulation ((7+) 
and the fcc-like modulation (q~) in the sense that when all are equally populated the cubic symmetry is preserved. 
Finally, at a temperature T^^, a martensitic-like transition to a ferromagnetic martensitic (tetragonal) phase (FMT) 
(e 7^ 0, g = g+, m 7^ 0, m,o 7^ 0) occurs. 

In figures g and Q we show the phase diagram as function of A'*. The intermediate transition (IT) is of first-order 
and it is represented by the dashed line (T/) separating, below the critical point (black diamond), the regions with g+ 
and q^ . For p = 2 (Fig. ||), this boundary is a vertical straight line and therefore cannot be crossed by sweeping T* 
at constant AT*. We advance here that this is a mean- field artifact. The exact treatment by Monte Carlo simulation 
will show that this transition line is always bent, even for p = 2. The mean- field solution renders a real first-order IT 
only for p > 2, as it can be seen in Fig. for p = 4. As an example, in figure o we show the temperature behavior of 
the order parameters for p = A and three values of K* around the critical value K*; K* < K* (a), K* ~ K* (b) and 
K* > K* (c) . The insets shows an enlarged view of the magnetization behavior around Tj . 

In figure || we show the location of the critical point (T/^) for p = 2 (AT* — 0) with respect to the magnetic (T^) 
and the martensitic (T^) transitions as function of the coupling parameter U*. One observes that, as the strength 
of the coupling decreases, the critical point approaches the martensitic line in such a way that the IT disappears well 
before U* becomes zero. This is consistent with some experimental observations ||18| indicating that strong coupling 
is required for the IT to appear. The same behavior is found for p = 4, as shown in Fig. M. We note that, in this 
case, the critical points correspond to different values of K* (AT* == — 2T*ln|). Figure g shows the location of the 



critical points (thick line) and the region of first-order IT in the U*-K* plane. For very weak coupling, the IT does 
not exists. For p = 4, as the strength of the magnetoelastic coupling increases, the IT exists for a larger interval of 
values of K* above X* . 

The magnetic transition is of second order for all values of model parameters studied. The MT is found to be 
first-order only for low values of the coupling strength for which e and q order simultaneously. Notice that this does 
not contradict the results obtained previously for the DBEG with U* — Q and K* > 0. Actually, for p ~ 4, it was 
found that the MT transition is discontinuous only for K* > ||l9[. Thus, the mean-field solution of the present 
model does not reproduce the evident first-order character of the MT in Ni2MnGa. At this point, this could be 
attributed to the insufficient accuracy of the mean- field treatment but, in the next section, we will see that Monte 
Carlo studies render the same feature. Thus, it is more plausible that it is due to the incompleteness of the model, 
addressed preferently to the study of the IT and related effects. Actually, our model description is done in terms of 
a single modulation strain whereas it is known |Q that the MT and the IT phases have different modulations. 

In conclusion, the main effect of the magnetoelastic coupling parameter U* is to generate a critical point (T/^, 
K* = — 2T* In |), between the magnetic (T^) and the martensitic (T^^) transitions. It is the end point of a first-order 
transition line Tj(U*) that, emerging from the martensitic phase boundary, separates two average cubic (e = 0) 
ferromagnetic (to ^ 0) phases with different modulation amplitude {fj = qo — q)'- q^ {fj < 0) and q~ (j) > 0). Provided 
the coupling strength is large enough, the IT exists for a limited range of K* (below the critical value) which, in turn, 
depends on U* . 

Before ending this section, we would like to show that the magnetoelastic coupling behind the model under discussion 
is consistent with the one considered in the Landau expansion (4), which in turn has been inspired by the experiments. 
Let us consider the simplest case of p = 2 (and K* = 0). Furthermore, we shall assume that to-o may be approximated 
by Too — (1 — q)m. Although this decoupling (see eq.(17)) is exact only for T* > Tj , it provides the first coupHng 
term between to and q. Indeed, from eq. (21) we see that the magnetoelastic contribution to the internal energy 
becomes of the form: 

-2C/*too(to - nio) ~ ~2U*q{l - q)rr? = ^^m^ + 2U*m^f, (26) 

where fj = qo — q=^ — q. As it was already discussed at the end of section (III), the first contribution represents a 
simple correction to the Curie temperature, and the second is the magnetoelastic coupling. 

V. MONTE CARLO SIMULATIONS 

In this section we solve numerically the model (^2|) by using Monte Carlo simulation techniques |Q. Our objec- 
tive is to find system configurations ({ci} , {•S'i}) distributed according to the canonical ensemble probability. The 
corresponding equilibrium simulations have been carried out using the standard Metropolis algorithm. The changes 
in the ai and Si variables are proposed independently and accepted or rejected according to the single-site transition 
probability W = min{l,e~^^ 1^ }. We have used a 2d-square lattice with N (= L^) sites subjected to periodic 
boundary conditions. Different lattice sizes ranging from L = 20 to i = 100 have been studied. The unit of time is 
the Monte Carlo step (MCS) and consists in N attempts of changing the Ui and Si variables. The simulations have 
been carried up to '--^ 30 • 10'^ MCS per site. Runs have been performed starting from two initial conditions: (i) a 
perfect FMT phase (cr^ = 1,5"^ = 1, i = 1, 2, ••• TV) and (ii) a perfect FC phase (oi = 0,5"^ = 1, i = 1, 2, • • • A^). This is 
very convenient in order to detect metastability and hysteresis when crossing first-order transition lines. Notice that, 
from the mean-field solution presented in the previous section, we already have an idea of the range of the space of 
parameters we have to explore. Accordingly, we shall fix J^ = 4.0 and use different values of V* < 0. Concerning the 
degeneracy factor we restrict the Monte Carlo simulations to p = 2. Nevertheless, we have verified that other values 
of p > 2 render qualitatively similar results. Most of the simulations have been performed at fixed values of the model 
parameters {U* and K*) and sweeping the temperature T*, but few have been performed at fixed T* and sweeping 
the parameter K*. 

The different quantities measured after each MCS are: the internal energy 7i*, and the order parameters to,, e and 
q. These quantities have been averaged over ~ 200 configurations taken every 100 MCS and discarding the first lO'' 
MCS for equilibration. Such averages will be denoted by (• • •). We have computed (7i*), {\m\),{q),{\e\). Moreover, the 
specific heat and the susceptibilities associated with the fiuctuations of the order parameters have also been measured: 

N ) T* 



c*-[^A "' '^.r' (27) 



_ (to^) - (|to|)^ 

x™ Y* '-28J 



..^^^V^ (30) 

All these definitions correspond to intensive quantities. In many cases the specific heat c* has also been obtained 
from the numerical derivative [1 / N)d{TL*) / dT . The agreement between this and the estimation obtained from eq. 
(27) gives confidence that the equilibration times used are appropriate. From the behavior of the order parameters, 
the specific heat and the corresponding susceptibilities the phase diagram can be obtained. The phase transitions 
associated with e, q and m have been determined from the location of the peaks in either the specific heat or in the 
corresponding susceptibilities p3[ . This method is more accurate than to look for singularities directly on the behavior 
of the order parameters. Moreover, we have checked whether of not the peaks correspond to a true phase transition 
by studying their dependence with increasing the system size L. As an example, in Fig. we show the temperature 
dependence of the specific heat for U* — —3.5, K* — 0.15 and four different system sizes {L — 10,25,50,100). The 
smooth peak at T* ~ 4.3 does not correspond to a true phase transition since it does not exhibit scaling behavior. 
The second order phase transitions (at T^j ~ 2.25 and T^ ~ 5.2) exhibit peaks which shift and become narrower 
and higher as one increases the system size L p3]. Besides, first-order phase transitions (Tf '--' 3.9 ) exhibit sharp 
discontinuities which, although they can also increase, neither shift nor become smoother with increasing L. In this 
last case if L is not very large or averages are not taken for long enough MCS, hysteresis may appear. 

Figure ^^ shows three sections of the phase diagram as a function of K* corresponding to three different values of 
the coupling parameter: U* — 0.0(a), U* = —2.5(b) and U* = —3.5(c). We notice that, for each value of U* , the 
magnetic (T"^) and martensitic (T^) transitions are almost independent oi K* . The overall conclusion emerging from 
Fig. |lO| is that the premartensitic effects are more important the larger the strength of the magnetoelastic coupling 
U* is. As it was anticipated by the mean-field calculations, its main effect is the showing up of a critical point (black 
diamond), and a first-order transition line (dashed line) separating two FC phases, g"*" and q~ , with different values of 
the modulation amplitude f]. Contrarily to the mean-field solution, for p = 2 , the Tj[K*) is now bent due to {(JiOj) 
fluctuations, as will be discussed in the next section. 

The second interesting point manifested by the Monte Carlo results is the existence of large fluctuations close to 
the critical point which, we stress, do not correspond to true phase transitions. These are revealed by anomalies 
in the response functions defined in equations (27-30). In the case of the specific heat, such anomalies appear in 
the form of smooth peaks that become difficult to resolve as we move away from the critical point. In figure [ll| , 
which is an enlarged view of figure |l0|(c), we denote the position of such smooth peaks by two dotted lines (with 
white squares) that, from the critical point, extend towards both sides of the FC phase. We have also indicated the 
metastability limits associated with the IT (points inside triangles). Actually, the position of the IT (denoted by 
black points along the dashed line representing Tj{K*)) is determined as the middle point of these limiting lines. 
The metastability limit points have been obtained by performing some runs at constant T* (< Tj*,) and sweeping K* 
(either increasing K* from the g"*" phase or decreasing K* starting from the q~ phase) and some others at constant 
K* and increasing the temperature from the MT phase. The first-order IT, as will be discussed below, is not found 
by decreasing temperature. 

Keeping Fig. nil in mind, now we shall study the temperature behavior of the specific heat at constant K* 
and U* = —3.5, above and below the critical point. The corresponding results are shown in Fig. ^ Different 
phenomenology may be observed when, increasing the temperature, we move from the bottom to the top of the figure. 

a) For K* — —0.1, a martensitic transition and a magnetic transition with no sign of an intermediate transition. 

b) For K* = 0.12, a martensitic transition, an anomaly (indicated by t) due to the proximity of the critical point 
and a magnetic transition. As we mentioned before, this anomaly is due to fluctuations and appears when 
crossing the dotted line (with squares) in Fig. |ll|. 

c) For K* — 0.14,0.15 and 0.18, an additional peak (to the left, also indicated by t) due to the intermediate 
transition shows up. As K* increases one observes that the fluctuation peak gets smoother (since we move away 
from the critical point) while the IT shifts towards lower temperatures. The entire temperature behavior for 
K* = 0.15 has been previously discussed in Fig. |^. 

d) Finally, for values oi K* even larger {K* — 0.25), the anomaly due to critical fluctuations has almost disappeared 
and only the peaks associated with the three phases transitions are clearly revealed: the martensitic, the 
intermediate (or premartensitic) and the magnetic transitions. 

It is difficult to be exactly at the critical point (that we have estimated K* ~ 0.130±0.005) but it would correspond 
to the value of K* at which the anomaly peak begins to split up. 
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It is interesting to note that such double peak behavior found around the IT has been observed experimentally in 
Ni2MnGa. Figure [l^ a shows an example corresponding to MC simulations with p = 2, U* = —3.5 and K* = 0.14 The 
two peaks are found for heating runs only. Due to the metastability of the first-order intermediate phase transition the 
transition to the q+ phase it is not found when cooling. A similar behaviour is found when performing calorimetric 
measurements as it is illustrated in figure n3b. The corresponding experimental details can be found in Ref. UM- The 
lines correspond to thermograms obtained by heating and cooling (as indicated by the leanning arrows). Actually, 
Fig. n3b corresponds to an enlargement of figure 1 in Ref. UW that has been reproduced with permission of the 
authors. Here, additional calorimetric runs are shown in order to reveal the systematic character of such double peak 
obtained when heating. Indeed, in the original published figure the two peaks are almost unobservable and were not 
considered by the authors who treat both peaks as a single one. When comparing Figures [l3| a and [l3p one observes 
that the highest peak occurs in different order. We do not have an explanation for this yet, but it could be related to 
the calorimeter inertia. 

Besides the behavior of the specific heat, it is also instructive to look at the behavior of the susceptibilities. Figure 
|lj shows the evolution of Xm, Xg ^nd Xe with temperature for K* = 0.15 and U* — —3.5 (this is one of the cases 
shown in Fig. |lj and discussed in point c)). The smooth peak in the specific heat (marked by t) is associated only 
with fluctuations of q, while at the IT (large peak) the discontinuity in Xq (modulating strain) is accompanied by an 
increase of the fluctuations Xe (tetragonal homogeneous distortion). 

Very recently, measurements of magnetic and transport properties of NiMnGa alloys have been reported [ p5[ . 
Magnetization measurements as a function of temperature, for small values of the applied field, reveal the existence 
of premartensitic anomalies at two different temperatures (separated '^ 20°). Only the low temperature anomaly is 
found when performing resistivity measurements. Nevertheless, it has been suggested pq] that both correspond to 
different true pre-martensitic transitions. In the light of the present results we suggest that the high temperature 
anomaly could be a signature of the critical fluctuations. We agree with the authors that more careful experimental 
studies are needed in order to clarify the results. 

VI. DISCUSSION AND CONCLUSIONS 

The phase diagram of the model exhibits, qualitatively, the same features in both mean-fleld theory and Monte 
Carlo simulations. A given value of the ferromagnetic interaction parameter J^ > determines the distance between 
the magnetic r,*^ and the martensitic T^j (< T*^) transitions. For appropriated values of the parameters K* and U*, 
the following phases are found, from high to low temperature: paramagnetic cubic (PC), ferromagnetic cubic (FC), 
ferromagnetic intermediate or premartensitic (FPMT) and ferromagnetic martensitic (FMT). The change from the 
FC to the FPMT phases occurs below a critical point (T/^, K*) and it takes place trough a true phase transition 
only for K* > K*. The existence of this intermediate phase depends on both K* and U* . First one needs the 
magnetoelastic coupling U* {—J^ < U* < 0) be strong enough. We obtain that Tj^ decreases with U* whereas Tj^ 
remains almost unaltered so that the critical point disappears, on the T^j line, well before U* reaches the value zero 
(Fig.^). Moreover, provided the U* is adequate, the IT exists for a limited range of values of K*{> K*) across the 
flrst-oder transition line Tj{K*). The corresponding order parameter is the modulating strain amplitude and changes 
from high to low values as one decreases the temperature across the IT. According to the theory of the harmonic 
thermal vibrations in a crystal, this is consistent with the behavior observed for the phonon frequency |14| . Moreover, 
in our results, the IT is accompanied by a jump in the magnetization, only visible for p > 2. Since this happens in 
both mean-field and Monte Carlo solutions we conclude that this has to do with the coupling rather than with the 



fiuctuations. Experimentally, a jump in the magnetization has not been detected so far [El 23 1. 

Some differences are obtained between the Monte Carlo simulation and the mean field solutions. First, in the 
simulations, the Tj{K*) line always bends towards increasing K* for p >2 whereas the mean-field solution gives, for 
p — 2 a. perfect vertical boundary at K* = 0. This is due to {ciCTj) fluctuations and may be understood as follows. 
The internal energy of a pure q~ phase {ui — 0, Si — 1) is given by (see eq. (12)): 

E- = -K*2N - J*,2N, (31) 

while the energy of a pure q'^ phase (ui — ±1, at random. Si = 1) is: 

^+ = - < ^ a,aj > -J^2N. (32) 

nn 

In the mean-fleld approximation the term between < .. > is neglected and therefore the condition E^ — E^ gives 
K* = 0. In Monte Carlo simulations, the fluctuations are present. We can estimate its value by taking the results 
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from the standard Ising model and therefore — < J2nn ^i'^i >— ^isingi"^*)- This function is zero only at T* = oo. 
For T* = takes the value —2N and increases monotonously. In the Ising model, at T* = T* ~ 2.27, Ej^^ = —V^. 
Therefore, the condition E^ = E+ gives a transition line at K*(T*) ~ -E}^^„g{T*)/2N. It bends towards K* > 
values when T* decreases. Remember that in the mean-field approximation a similar behavior is found only when the 
degeneracy factor of the 0-state is p > 2. 

Another point illustrated by the Monte Carlo simulations is the important role played by the fluctuations in 
describing premartensitic effects. In figure |lj we show, schematically, the K* — U* section of the phase diagram 
as it is obtained from the numerical results presented in the previous section. The triangle defines the region with 
a first-order intermediate (premartensitic) transition and it is limited by the line of critical points (thick line) and 
the line where the IT disappears on the MT line (thin line). The shadow region denotes the zone of large (critical) 
fluctuations and it has been estimated from the anomalies in the speciflc heat. It is interesting that, apart from the 
region of fluctuations, the mean-field solution renders a qualitative similar phase diagram (Fig.||) although shifted to 
negative values of K* . 

Experimentally, the interplay between U* and K* , required for the IT to occur, is determined by the composition 
of the sample. Indeed, premartensitic effects in the NiMnGa alloy have only been reported for compositions around 
the stoichiometry (Ni2MnGa). Among them, the most important and common to all the samples, is a significant 
softening of the ^ [1 10]TA2 acoustic mode with decreasing temperature, accompanying the formation of an intermediate 



structure that, while preserving the cubic symmetry, is transversely modulated [10 36| along the [110] direction with 
wavevector ^ . Moreover, this intermediate structure may appear through a true phase transition [|l^,|lj[ or not [ p7[ . 
We point out that both behaviors are compatible with the present results. In simple words, the samples showing the 
IT would fall inside the triangle in Fig. |l^ while the others not. More precisely, in the samples showing the IT the 
strength of the magnetoelastic coupling (f/*) is enough to freeze completely the anomalous phonon with degree of 
softness related to K* . For deeper discussions it is imperative to have a better understanding about what is behind 
the composition dependence giving rise to the different behaviour in the samples. In particular, experiments in order 
to compare the different degree of softness of the anomalous phonon are required. Furthermore, it is important to 
know the strength of the magnetoelastic coupling. For this, measurements in both the ferro and the paramagnetic 
phases of the temperature behaviour of wf are needed. An study of the characteristics of the kink around the Curie 
point would be very helpful. Finally, we point out that samples in the critical region (shadow region in Fig nS) should 
exhibit a significant increasing of diffuse scattering when decreasing the temperature below the Curie point. 

The present model renders, independently of the technique used to solve it, a MT which is continuous for the range 
of model parameters studied. This would be a serious setback in case we where interested in the properties of the MT 
itself. For this, models such as discussed in references [g[ and [g^ could be more appropriate. Here, we have developed 
a model with the aim to focus on the study of the IT and related premartensitic effects in Ni2MnGa. It is based on 
the assumption, sustained by the change of the slope -^ at the Curie point, that the magnetism causes the freezing of 
the incipiently unstable |[110] TA2 phonon and the splitting from the homogeneous strain. With this point of view, 
the intermediate phase is a precursor of the MT. On the other hand, in view of the fact that the actual modulation 
of the martensitic phase is different from that of the intermediate phase, some authors have claimed that both phase 
transitions have to be regarded as independent [|37| . If so, additional coupling terms between the modulation and the 
homogeneous strains are required in order to produce a change in the modulation at the martensitic transition. In 
this sense, it has been observed that the dip in the TA2 branch shifts under an external unaxial stress [Q. However, 
more along with the point of view adopted here, Sthur et al. [0[ have demonstrated that the [|-|0] characterizing 
the intermediate phase becomes a vector [0.38, 0.38, 0] in the tetragonal phase, very close to the [l|0] expected for a 
5-layered modulation. 

From the above discussion is follows that, in spite of the appealing results obtained, the present model needs to 
be improved in order to reproduce the whole scenario of the structural transition in Ni2MnGa. Finally, it is worth 
mentioning that, by increasing the strenght of the magnetoelastic coupling U* , it is possible to extend the present 



study to samples for which the MT takes place in a paramagnetic matrix |39 
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TABLE I. Identification of the different phases of the model and their corresponding abbreviated notation. The phase q 
and q~ axe both average cubic phases and undistinguishable above the critical point. For convenience, the q~^ phase is also 
called premartensitic FPMT. Note that qo = — r, takes the value 1/2 for p = 2. 



Name 


e 


m 


<? 


Paramagnetic cubic PG 








- 




pure cubic 





/o 


qo 


Ferromagnetic cubic EG 


average cubic (g~) 


< go 




Premartensitic FPMT (g+) 


> go 


Ferromagnetic martensitic FMT 


/o 


/o 


> go 



13 



p=2, U =-3.5, K =0.0 



(7) 
DC 
LU 

I- 
LU 

< 

DC 
< 

CL 

DC 
LU 
Q 
DC 
O 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



\" 


' 1 ' 


1 


\ 8 


+\ \ 

q \ \ 


- 


1 


/ \ 


" q 


" 


qy \ 






/ \ 


m 


1 


, 1 , 





0.0 0.5 1.0 1.5 2.0 2.5 



3.0 



FIG. 1. Mean-field temperature evolution of the order parameters for J* 
line, that in this case (p = 2) corresponds to K* = 0. 
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FIG. 2. Mean-field, temperature evolution of the order parameters for J^ = 4.0, U* 
line K*{T*) = -2T*ln{p/2). 
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FIG. 3. Section of the phase diagram obtained by mean-field techniques for J^ = 4.0, U* = —3.5 and p = 2. Continuous 
hnes stand for second-order phase transitions and dashed hues for first-order ones. The black diamond shows the location of 
the critical point. The labels indicate the different phases as explained in Table I. 
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FIG. 4. Section of the phase diagram obtained by mean-field calculations for Jj^ = 4.0, U* = —3.5 and p = 4. Continuous 
lines stand for second-order phase transitions while dashed lines for first-order ones. The black diamond shows the location of 
the critical point. The labels indicate the different phases as explained in Table I. 
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FIG. 5. Mean-field temperature evolution of the order parameters for J^ ~ 4.0, U* — —3.50, p = 4 and three diflferent values 
of K": (a) K* = -2.0 < K*, (b) K* = K* and (c) K* = -1.85 < K*. The insets show the detail of the behavior of the 
magnetization m at the intermediate transition. 



18 



4.0 



3.0 



2.0 



1.0 



0.0 



p=2, K =0 




FMT 



0.0 



-1.0 



-2.0 

* 

U 



-3.0 



-4.0 



FIG. 6. U*-T* section obtained from mean-field calculations for J^ = 4.0, p = 2 and K* = 0. Continuous lines stand 
for second-order phase transitions. The labels indicate the different phases as explained in Table I. The intermediate region 
indicated by FPMT is a phase separation region with coexistence of the two phases with q^ and q" . The thin dashed line 
indicates the section of the phase diagram shown previously in Fig. 3. 
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FIG. 8. U*-K* section of the phase diagram obtained from mean-field calculations for p — 2 and p — 4. The location of the 
critical points are shown with a thick continuous line. The thin line is estimated from the value of K* at which the IT and 
the MT meet. The region inside the triangle denotes the range of the model parameters for which there exists the first-order 
intermediate transition (IT) in the case of p = 4. 
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FIG. 11. Detailed section of the phase diagram (Fig. 10(c)) obtained by MC simulation with U* = —3.5 and Jm ~ 4, showing 
the details of the metastability regions and the fluctuations in the neighborhood of the critical point. Symbols and lines have 
the same meaning as in Fig. 10. Besides, we have indicated the position of the large fluctuation peaks with squares and the 
metastability limits of the first-order transition lines with points inside triangles. The orientation of the triangles indicates 
the direction of the MC simulation runs performed in order to locate each metastability limit. Note that we have performed 
increasing T* runs (A) and increasing ( >) and decreasing (<|) K* runs. Dotted lines are guides to the eyes 
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FIG. 12. Specific heat evolution, corresponding to U* — —3.5 and different values of K* as indicated. Tfie arrows on eacfi 
plot indicate the peaks corresponding to the IT and the large fluctuations due to the closeness of the critical point. 
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FIG. 13. Details ol the specific heat evolution at the premartensitic phase transition showing the double peak effect and hys- 
teresis, (a) Results from MC simulation of the present model corresponding to U* — —3.5 andA'* = 0.14 and (b) experimental 
data extracted from Ref. 11 
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FIG. 14. Comparison of the specific heat behavior (a) and the susceptibilities (b). In (a), the two curves correspond to the 
estimation of c* from the derivative of the energy (•) and from the fluctuations as defined in the text (o). In (b) the three 
different curves correspond to the fluctuations of the three order parameters e, q and m, as indicated by the legend. 
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FIG. 15. Space parameters U* —K* obtained with MC simulations, indicating the regions with IT and with large fluctuations. 
The thick line indicates points with true critical behavior. 
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